Spin freezing by Anderson localization in one-dimensional semiconductors 
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One-dimensional quantum wires are considered as prospective elements for spin transport and 
manipulation in spintronics. We study spin dynamics in semiconductor GaAs-like nanowires with 
disorder and spin-orbit interaction by using a rotation in the spin subspace gauging away the spin- 
orbit field. At a strong enough disorder spin density, after a relatively fast relaxation stage, reaches 
a plateau, which remains a constant for long time. This effect is a manifestation of the Anderson 
localization and depends in a universal way on the disorder and the spin-orbit coupling strength. 
As a result, at a given disorder, semiconductor nanowires can permit a long-term spin polarization 
tunable with the spin-orbit interactions. 

PACS numbers: 72.25.Rb,72.70.+m,78.47.-p 



I. INTRODUCTION 

The main idea of spintronics - the design and applica- 
tion of devices controlling not only the charge dynam- 
ics but also the electron spin evolution - can be use- 
ful for information storage, transfer, and manipulation 
technologies^— Possible realizations of spintronics de- 
vices can be based on semiconductor nanowires^— for 
quasi-ballistic electron transport, coherent transmission 
of information, and spin control. These systems attract 
a great deal of attention due to a clear interplay of trans- 
port and spin-orbit (SO) coupling physicsJ^ — 

This control faces the problem of inevitable spin relax- 
ation due to the coupling of electron spin to environment 
through SO coupling. As a result, the factors determin- 
ing the spin relaxation rate become of crucial importance. 
Two limiting cases of spin relaxation are well understood. 
For the itinerant electrons spin relaxation in mainly de- 
termined by the Dyakonov-Perel' mechanism, that is by 
random precession of electron spin due to the random in 
time electron momentum. 

A different approach should be applied for elec- 
trons localized in a regular external potential form- 
ing quantum dots promising for quantum information 
applications^ Here momentum is not a well-defined 
quantity, and the momentum-dependent splitting re- 
quired for the Dyakonov-Perel mechanism vanishes. As 
a result, spin relaxation through SO interaction requires 
phonon-induced coupling of different orbital states of the 
localized electron and nonzero external magnetic fields 
In the absence of magnetic field and spin-orbit coupling, 
spin relaxation can occur due to the hyperfine coupling of 
electron spin to spins of lattice nuclei In both cases, the 
initial spin polarization goes asymptotically to zero. The 
characteristic timescale of spin relaxation of electrons lo- 
calized in quantum dots is expected to be several orders 
of magnitude longer than that of itinerant electrons. 

While these two limits of free and strongly localized 
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FIG. 1: (Color online) Semiconductor nanowire with random 
impurities shown as filled circles. Although we consider a 
one-dimensional electron motion, impurities can be randomly 
distributed over the cross-section of the wire. 



electrons are well understood, the interplay of disorder- 
induced localization and spin relaxation of itinerant elec- 
trons remains an open question although some aspects 
of the problem have been addressed^ — The questions 
here are (i) how the localization forms the spin relax- 
ation, and (ii) whether the initially prepared spin density 
relaxes to zero. As a nontrivial example of this interplay 
we mention that weak localization of two-dimensional 
electrons leads to a long power-like rather than exponen- 
tial spin relaxation ! 24 i 25 Here we analyze this problem 
for the one-dimensional system, providing, on one hand, 
the basic example of localization physics in a random 
potential ^ 26 ! 27 and, on the other hand, an example of a 
system, where spin-orbit coupling can be gauged away 
by a SU(2) transformation. 

This paper is organized as follows. In Sec. II, we show 
how to treat spin relaxation in one-dimensional systems 
with the gauge transformation and introduce the tight- 
binding Hamiltonian for the model. The spin dynamics 
will be analyzed by a numerically exact calculation in 
Sec. Ill, where we show that spin density does not re- 
lax to zero, in contrast to what expected. In addition, in 
this Sec. Ill we study how asymptotic value of spin polar- 
ization depends on the disorder and spin-orbit coupling. 
Conclusions summarize the results in Sec. IV. 
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II. MODEL 
A. Hamiltonian and gauge transformation 

The investigated structure is a quantum wire extended 
along the x axis, as shown in Fig. [TJ The total Hamilto- 
nian has the form 



H 



— (k x A x ) 
2m 



U(x) 



ma 



(1) 



where A x = —maa y /h stands for the Rashba coupling 28 
with the strength a, a y is the Pauli matrix, k x is the elec- 
tron wave vector, and m is the effective mass. The Dres- 
selhaus coupling 2 ^ is obtained with A x = —m(3a x /h 2 , 
where f3 is the coupling constant. Without loss of gener- 
ality, we concentrate here on the Rashba coupling, which 
can be changed on demand by applying external electric 
field across the structured 

The SO interaction can be removed from H in Eq.(pQ) 
through a gauge transformatio n 31 ! 32 with a SU(2) spin 
rotation: S = exp (— ixa y /2^) , where £ = h 2 /2ma is 
the spin-precession length. After this transformation the 

system Hamiltonian has the form: H = Ti 2 k 2 /2m J r U(x). 

Since for the Hamiltonian (pQ), <j y is the integral of mo- 
tion, the spin density component along the ?/-axis is time 
independent. A nontrivial dynamics of the transformed 
spin occurs for the 7 = (x, z) spin components (s 7 (x,t)) 
and can be expressed in terms of the spin diffusion 



(s 7 (x,t)) 



J D^P(x- 



x',t)(sp(x',0))dx', (2) 



where D 1 ^(x 1 t) is the exact disorder-dependent one- 
dimensional spin diffusion Green's function. In a non- 
magnetic system without SO coupling D 1 ^(x 1 t) = 
8 1 pD(x, t) is diagonal in the spin subspace. As a result of 
the gauge transformation, the uniform density dynamics 
is determined by only the Fourier component 25 



D(q,t) 



dxe- iqx D(x,t) 



(3) 



with q = 1 /2£ and Eq. (J2]) simplifies for the physical mea- 
surable spins as (s 7 (£)) = (s 7 (0))Z)(l/2£, t). Here we 
will use a similar, however, somewhat different approach 
based on numerically exact analysis of the direct time 
evolution of the initial spin-polarized states. It will be 
shown that the resulting spin dynamics has unexpected 
features, including a long-time plateau in the spin polar- 
ization. 

The eigenfunctions of H can be taken in the form 
ijj{x) = ifj(x) |1) and ip(x) = ip(x) | — 1), where |±1) are 
the eigenstates of <j z with the corresponding eigenvalues. 
The eigenstates of H, <j>{x) can be obtained by spin ro- 
tation of the ip(x) \a) states. For example, with spin- up 
initial state ip(x) |1) one obtains: 



(x) = i/j(x) 



(4) 
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FIG. 2: (Color online) Site dependent components of 
(j) N / 4 (xn) for a qualitative description of entanglement in- 
duced by the gauge transformation for (a) a — 0, (b) a = 
0.125 x 10" 6 meVcm, and (c) a = 10" 6 meVcm (U = 55 
meV). The solid and dashed lines represent |1) and | — 1) com- 
ponents, respectively. 



The spin dynamics and spin relaxation in the system, as 
it will be shown below, is solely due to the entanglement 
of spin and coordinate in Eq.(J2]). 



B. Tight-binding model and disorder 

We perform numerical analysis using the tight-binding 
model, employing the approach similar to Refs.[16.33]. 
The one-dimensional electron gas is sampled with TV = 
2 13 (8192) grid points x n = nl, where 1 < n < N and 
I is the effective lattice constant with periodic boundary 
conditions^ The effective hopping matrix element be- 
tween two nearest neighbors is chosen as t = 50 meV, 
and the kinetic energy is E(k x ) = 2t(l — cos(k x l)). As a 



result, the eigenenergies span the range of [0, 200] me Vs. 
The distance between two neighbor grid points becomes 
I = h/\/2mt = 3.37 nm to satisfy the electron effective 
mass m = 0.067 mo in GaAs semiconductor with mo 
being the free electron mass. 

The random potential U n = U(x n ) uniformly spans 
the range [— Uq/2, Uo/2] with the white noise correla- 
tor (U(x ni )U(x n2 )) = (U 2 )5 nun2 , where (U 2 ) = U 2 /12. 
The effects of disorder can be approximately charac- 
terized through the energy-dependent momentum relax- 
ation time te, which we define as Ti/te = (U 2 )Ive-, where 
ve = \fm/7rhV2E is the density of states per spin com- 
ponent. The_£esulting mean free path £e = vete-, where 
ve = \/2E/m is the electron speed and the correspond- 
ing diffusion coefficient De = v e te- 

In this representation the eigenstates of H and H form 
basis sets, {V^} and {<^} respectively, where 1 < i < 2N. 
For the same z, these two sets are related by the local spin 
rotation S. We assume that = ipi(x n ) |1) for 1 < i < N 
and i) i i/>i-N(x n ) for N < i < 2N. 



III. SPIN DYNAMICS 

We study dynamics of initial r i\) i states with 1 < i < TV, 
corresponding to the evolution upon instant switching of 
the SO coupling. The time dependence can be expressed 
with the spectral decomposition as: 



= 1,2N 



-itei /h 



(5) 



/ gQ gQ \ 

(cr z (t))j = (ipj (t) cr z ipj (t)) is determined by the spec- 



where = (4>i\ipj), and e% are the corresponding 
eigenenergies. The spin component expectation value 

— so, 

trum and eigenstates of the system. 

In order to give an idea of the entanglement induced 
by SO coupling, we present in Fig. [2] the evolution of 
4>N/A( x n) state with the increase in the spin-orbit cou- 
pling. At a = 0, we obtain a product state 4>N/A( x n) — 
ipN/4:( x n) and with the increase in a entangled states 
are formed. The overlap of 4>i(x n ) and ipj(x n ) eigenstates 
is characterized by two sets of matrix elements a^-; for 
example, for 1 < j < N : 

dij = y^CQS (j^J ^i( x n)^j( x n), l<i<N (6) 

dij = 5Z Sin ( ) ^i( x n)^j(Xn), N <i <2N. 
n \ s / 

The behavior of presented Fig. [3] demonstrates that 
for given j it has nonnegligible values only in a certain, 
rather narrow, range of i. 

To illustrate the role of the random potential, we con- 
sider as examples weak (Uq = 15 meV, Uq <C t) and 
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FIG. 3: (Color online) Absolute values of ciij around the ini- 
tial spin-up state ip N / 4 ; here a = 10 -6 meVcm (strong SO 
coupling) and Uo=55 meV (strong disorder). 
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FIG. 4: (Color online) Inverse participation ratio £ for the low 
part of the energy spectrum; gray (red) circles denote strong 
disorder (£/ =55 meV) and black circles denote weak disorder 
(Uo=15 meV). Since even for Uo=55 meV we obtain (< 1, 
the localized states are distributed over many lattice sites, 
confirming applicability of the tight-binding Hamiltonian for 
the localization problem. 



strong (Uq = 55 meV, Uq > t) disorder. For free elec- 
trons in state j = 7V/4 and E = 31 meV, the result- 
ing Ti/te is about 0.1 and 1 meV, respectively. For a 
free electron with the energy E w 20 meV the velocity 
ve ~ 3.5 x 10 7 cm/s, the mean free path £e ~ 2.5 x 10 -5 
cm (Ti/te = 1 meV), and the corresponding diffusion co- 
efficient De ~ 10 3 cm 2 /s. These parameters provide an 
effective integral characteristic of the disorder and corre- 
spond to realistic parameters of the wires, which, how- 
ever, can strongly vary from sample to sample and from 
experiment to experiment. 

The effect of localization by disorder is seen in the in- 
verse participation ratio 3 -^ (IPR) d = J2 n I^O^n)!- The 




FIG. 5: (Color online) Time-dependent polarization in the 
weak-disorder regime (U = 15 meV). The initial bins are cen- 
tered at the states (a) N/4 (bin width 6.8 meV), (b) N/8 
(bin width 3.7 meV), and (c) and N/16 (bin width 2.1 meV) 
with energies decreasing in the same order. The curves 
for SO couplings 0.125 x 10" 6 meVcm, 0.5 x 10" 6 meVcm, 
and 2 x 10 -6 meVcm are drawn with circles, triangles, and 
squares, respectively. Note that after the relaxation stage the 
spin density remains a finite constant. 



IPR calculated for the low-energy spectrum is presented 
in Fig. |4j As expected, the degree of localization increases 
with Uo and this effect is more pronounced for the elec- 
trons with lowest energies. In contrast to the results of 
Ref.(23[, the IPR in this system does not depend on the 
SO coupling. We now study the effects of disorder and 
spin-orbit coupling on the average spin dynamics of a bin 
of 256 initial spin-up states and 8 realizations of the ran- 
dom potential. The statistical error of this approach is, 
therefore 1/^2048=2.2%, making the results statistically 
representative. 

We take three example bins with three different de- 
grees of localization. The bins are centered around the 
spin- up states ^n/& an d ^n/iqi whose IPR val- 

ues increase in the same order (energies decreasing, see 
Fig. 3]). The calculated bin- and potential realization- 
averaged spin dynamics is shown in Figs. [5] and [6j re- 
vealing strong influence of the disorder-induced spatial 
localization of states. Physically, collisions of electrons 
with impurities force electron spin to frequently reverse 
the precession direction. In the classical picture, this 
leads to a long Dyakonov-Perel' spin relaxation. If the 



FIG. 6: (Color online) Time-dependent spin polarization in 
the strong-disorder regime (?7o=55 meV) with the same no- 
tations as in Fig. El (a) N/4 (bin width 7.2 meV), (b) N/8 
(bin width 4.5 meV), and (c) JV/16 (bin width 3.7 meV). Note 
that for a = 0.125 x 10" 6 meVcm the spin is almost constant 
in time, thus suitable for spin-based operations. 



quantum effects of localization are important, the result- 
ing effect is the "freezing" of the electronic spin. As one 
can see in Figs. [5] and [6l the electron spin density re- 
laxes for ~ 5 ps and then remains constant in time for 
infinitely long (beyond 0.2 ns in our computation). As 
expected, the spin polarization plateau is higher (i) for 
localized states and (ii) for weak SO interaction. Al- 
most time- independent spin states are achieved e.g., at 
Uo = 55 meV and a = 0.125 x 10 -6 meVcm. 

To gain insight into the problem, we study the depen- 
dence of asymptotic spin density on SO coupling and the 
localization of electrons in more detail. The long-term 
densities are plotted in Fig. [7] against parameter £(£)■ 
This parameter combines the two factors determining 
the spin dynamics, SO coupling and spatial localization, 
where (() is averaged over 256 bin states and 8 realiza- 
tion of the random potential. The given values follow a 
universal dependence indicating a unique trend for long- 
term spins against SO coupling and localization through 
disorder. This trend corresponds to a fast increase in the 
asymptotic steady polarization for < 1 and a smooth 
increase and saturation for > 1. These results can be 
understood as follows. To show an efficient spin dynam- 
ics, the electron should move the distance of the order 
of 7r£. Therefore, the spatial spread of the correspond- 
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FIG. 7: (Color online) Long-term relative polarization as a 
function of for three different degrees of localization. Pa- 
rameter £ is modified by changing the coupling constant a. 



ing states should be larger than 7r£. With a stronger 
localization, the spread and the overlap decrease leading 
to the universal behavior shown in FigJTJ Qualitatively, 
in the "clean" <C 1 regime the spin relaxation has 
the Dyakonov-Perel' mechanism either purely exponen- 
tial for VLeTe < 1 or a combination of oscillations and 
exponential decay if SIe^e > 1, where the spin preces- 
sion rate Qe — 2a\/2mE /% corresponds to the electron 
momentum at given energy E. 



culations show that a prepared spin density relaxes un- 
til reaching a plateau, directly related to the disorder 
and strength of SO interaction. In contrast to the ex- 
pected decay to zero, a long-time constant polarization 
plateau survives to infinite time. The asymptotic spin 
density has a universal dependence on the product of 
the inverse participation ratio and the spin precession 
length. In the absence of magnetic field, the hyperfine 
coupling to the spins of nuclei will lead to spin relaxation 
on timescales at least two orders of magnitude longer 
than the timescale of the plateau formation of the or- 
der of 10 psJ^ As the experiments on spin transport did 
not reveal electron-electron interaction effects^ here we 
have neglected them. Furthermore, whether there exists 
a range of parameters where the Coulomb forces can be 
strong enough to modify our results for localized states, 
remains to be investigated. 

An immediate consequence of this result is the ability, 
by choosing the desired Rashba SO parameter for a given 
wire, to produce and destroy steady spin states, which 
are of interest for spin-based operations. These results 
suggest that semiconductor nanowires can be used for 
coherent transmission and storage of information, manip- 
ulated by spatially and temporally modulated spin-orbit 
coupling. 
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